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Error estimates are given for ||i<j — 9?||, where Ui is a normalized eigenvector of a self- 
adjoint operator A and <p is a normalized approximating vector. The estimates contain 
upper and lower bounds to certain of the eigenvalues of A together with the terms {A<pj<i^ 
and (A<p,A(p). 

1. Introduction 

In this paper we give some estimates for the error in norm between eigenvectors of a 
self-adjoint operator and approximating vectors. Our estimates involve upper and lower 
bounds to certain of the eigenvalues togetlier with terms in the approximating vector itself. 
Estimates using upper and lower bounds alone have been given by Lowdin and Sliull ^ and 
Weinberger.^ Other estimates involving the same quantities used here have been given by 
Kato,^ Kryloff/ Trefftz,^ and Weinberger.^ Such error bounds are useful in the estimation 
of quantum mechanical expectation values, error estimates for approximate solutions of 
Au=f, and other applications. 

2. Simple Eigenvalues 

We suppose ^ to be a self-adjoint operator with domain 3)^ in a separable Hilbert space 
§ having inner product (UjV). A is assumed to be bounded below and to have at least the 
initial part of its spectrum consisting of eigenvalues of finite multiplicity. These eigenvalues 
are considered to be ordered in a nondecreasing sequence, Xi<X2<X3< . . ., in which each 
eigenvalue appears according to its multiplicity. We denote the corresponding orthonormal 
eigenvectors by Ui, U2, u^, . , . , 

We wish to estimate the difference in norm between an eigenfunction Ut corresponding 
to a nondegenerate eigenvalue X^ and a given vector (p in J)^. For convenience we wiU con- 
sider (p to be normalized but undetermined to the extent of a scalar multiplier of magnitude 
one. The norm difference between Ut and (p is given by 

\\(p—Ui\\ = i<p—Ui,(p—u^y^=^J2[l — Ee(ip,u^)]K 

If we suppose that the unspecified scalar multiplier of <p has been chosen so that ((p,Ui) 
is real, positive, and equal to \((p,Ui)\, then 

so that our problem is equivalent to the determination of a positive lower bound to the 
magnitude of the projection of cp on Ui. 
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In order to obtain our estimates we require upper and lower bounds to a finite number of 
eigenvalues. We shall designate a lower bound to the jih eigenvalue by Xj and an upper bound 
by X^; thus 

\'j<\j<\^. (1) 

A condition sufficient to guarantee that Xj be nondegenerate is that X"_i<^Xi and X^+i^X". 
Upper bounds are usually found by the Rayleigh-Ritz procedure, while lower bounds can 
often be obtained by a variety of procedures. 

Let us choose n>i. By the spectral theorem we have ^ 

{Acp,cp)= r Xd{Ex<p,cp) > f: \j\(<p,uj)\' +Xn+i[l- f: \{<P,Uj)\'l (2) 

J - 00 i = 1 j=i 

Smce the coefiicients of Xi, X2, . . ., X^+i are nonnegative, we may replace these eigenvalues by 
the lower bounds Xj, X2, . . ., X^+j to obtain 

{A<p,^) > ± m<p,u,)n\u^ fi- s i(^,%)i^1 • (3) 

i=i L ;=i J 

The inequality (3) may be solved for \((p,Ui)\^ in the form 

where S' means that the term for which j equals i has been omitted. This lower bound for 
\((p,Ui)\^ still involves unknown eigenvectors of A. This difficulty will now be overcome 
through use of the Schwarz inequality. For the arbitrary real number fXj we have 

[h-H?\ i^J,<P) P= I {[A-ixj]uj,,p) p= ! {uj,[A-,.M'< 1 1 [A-nM IS (5) 

SO that for fij 7^ Xy we have 

|(„„)p<I^=i^'. (6) 

We assume that ip is not one of the eigenvectors %i^ u^, . . ., i/y-i, 'Wy+i, . . .^ Un and minimize 
the right hand side of (6) with respect to ju;. For the minimizing value given by 



{Aip,Aip)-\{^A^,<p) 



the inequality (6) becomes 



\U„\\i<i {Aip,A^)-{Aip,ipy 



Let us define numbers 7^ by 

' \\-{Aip,ip), {A^,ip)<\\ 



7i= 
so that 



, X5<(^^,^)<X,- \^ (9) 

I, {A^,<p)-^, (^^,^)>X," 

7?<[(^W)-Xi?, i=l,2, ...,7i; OVi). 
Using these mequalities, the inequalities (8) can be replaced by the weaker ones 



7 In the case that A has only n eigenvalues before the first limit point X* of its spectrum, we replace X„+i by X* in (2) and the succeeding 
equations. 



Our lower Ixniiid for \{<p,Ui)\^ is obtained by using the inequalities (10) in (4); wo have * 

'^^''*--^l - \J,+,-X| XU-XI 'j^ (^M^)-(^^,^)^+7?" ^''^ 

When the right hand side of (11) is positive, it can be used to give an upper bound to ||v?— ^i||. 

3. Some Special Estimates 

In those cases when there are an infinite number of eigenvahies converging to a finite 
limit point or diverging to infinity it is possible to make other estimates of \{(p,Ui)\^ which 
contain infinite sums. When these sums converge and can be esthnated numerically the 
resulting inequalities may be useful in applications. We briefly sketch two such results. 

Let us consider the case when the initial spectrum of A consists of an infinite number of 
eigenvalues Xi<X2< . . . that converge to the finite limit point X*. From the spectral theorem 
we have 

(A^,^) > S XJ \(<p,uj)\'+xdl--j: I {<p,uM''\^ (12) 

d=l L j=i J 

wliere x] and X^ are lower bounds to X^- and X* respectively. If Xj is a simple eigenvalue and 
if X^>Xj, then (10) and (12) give 

\f ..\\2^ ^*-(^<Py^) {A<p,Aip)-{Aip,(py ^, \l-\] , 

\{<P,u,)\ > ^,_^, - ^j_^, Z. (A<p,A^)-{A<p,<py+yf ^^^^ 

Since the y/s are again given b}^ (9), the inequality (13) requires the knowledge of upper bounds 
to each X^ for which (A(p,(p) >\]. 

When A has a pure point spectrum diverging to infinity various estunates can be given. 
One of them starts from tlie Parseval equation 

7 = 1 

and leads directly to the estimate 

4. Degenerate Eigenvalues 

We now suppose that Xj is either an eigenvalue of multiplicity m, or one of m eigenvalues 
which lie close together. Let p and q be the largest nonnegative integers such that 

Xj \ Xj_p 
and 

X"^Xi+j; 



8 The inequality (U) has the following immediate generalization. If /(X) is a real valued function which is monotonically increasing and 
fdA) is the corresponding self-adjoint transformation, then 



\i<P> 



^ ' - fiK+i)-f0^i) /(^i+i)-/(^!) ill {A<p,A<p)-(A<p,<p)^-\-yf 



then m=p-\-q+l. One appropriate question that may be asked is ^'How closely is some 
eigenvector in the span of Ui-j„Ui-.p^i, . . ,,Ui+q approximated by a unit vector <pV' This is 
equivalent to finding a lower bound to the projection of (p on the span of ^^^-p, . . .,Ui+q, It 
follows in the same way as (11) that 

where n>i+q and 2' means that the terms for which j=i—p,i—p+lj . . .,i+q have been 
omitted. The estimates (13) and (14) can be extended in the same way. 
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